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We study force-free configurations of Abrikosov flux lines in the line-liquid and line- 
crystal limit, near the melting transition at H m . We show that the condition for zero 
force configurations can be solved by appealing to the structure of chiral liquid crystalline 
phases. 
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In high temperature superconductors, the extremely large ratio of the London pen- 
etration depth A to the coherence length £ suggests that the most important degrees of 
freedom are Abrikosov vortex excitations. The configuration of flux lines in applied cur- 
rents and fields thus becomes of great interest. It is therefore useful to construct a theory 
of the flux-lines themselves which may be used to study their conformations. 

We start by considering the London equation for a superconductor, which relates the 
current density j to the magnetic field B: 

Vxj^-j^B (1) 

where A is the London penetration depth. This is, of course, supplemented by Maxwell's 
equation: 

47T 

V x B = — j. (2) 

c 

These two equations predict much of the phenomenology of superconductors. In particular, 
if B is along the 5-axis, then the London equation predicts that screening currents will 
circulate in the xy-plane. If the magnetic field is confined within flux-tubes parallel to the 
z-axis then the currents will circulate about these confined regions. These flux tubes will 
form the Abrikosov flux-line lattice. 

Under an applied current, flux-lines will adopt a steady state configuration in which 
there is no net force. In the absence of pinning the Lorentz force per unit length on the 
vortices, F/£ = B x j, will be balanced by repulsive vortex-vortex interactions. In a 
"force- free" configuration [|l| it is necessary that j be parallel to B: 

a(x)B(x) = ^j(x) = V x B(x) (3) 

where the last equality follows from Maxwell's equation and ct(x) is a spatially-varying 
scalar. In a superconductor, the magnetic field is confined to be near flux tubes and along 
their tangents. If the flux-lines trace out the curves Ri(s) where s is their arclength, then 

m(x) = J ^^(R^-x) (4) 

i 

is the local tangent density of flux-lines. Since flux-lines cannot end we have 



V-m = 0. 
1 



(5) 



Using standard techniques for treating topological defects we have: 



[1 - A 2 V 2 ] B = $ m (6) 

where $o is the flux quantum. It is useful to decompose m as the product of a unit vector 
n and an areal density p, m = pn [[| . If a varies on a lengthscale long compared with the 
penetration depth then by applying the operator [l — A 2 V 2 ] to (|3]) we find: 

Vp x n + pV xiiK a ph. (7) 

If we consider the system near H m where the flux-lines are dense, we can take p ~ po, a 
constant. In this case (|5|) and (^) become: 

Vn = 

( 8 ) 

V x n = aii 

These two equations are familiar in the field of liquid crystals: the flux-lines will attempt 
to adopt a configuration with no splay (V-n = 0), no bend (n x [V x n] =0), but with 
twist n-V x n = a. We will pursue this analogy with liquid crystals. Of course, the 
flux-line density does not need to be uniform. The liquid crystal analogy will allow us, 
however, to consider a class of paradigmatic vortex configurations which do not require 
density variations and are thus of low energy. 

For simplicity, we consider a superconductor in a magnetic field, applied along the 
z-axis. The Abrikosov flux lattice can be modeled as an elastic medium [|J: 

^Lattice = ^ J d 3 x {cnu| + 2c 66 [u^ - u%] + c 44 (<9 2 w) 2 } , (9) 

where u is the two-dimensional displacement vector (perpendicular to the average flux line 
direction), Uij is the two-dimensional strain tensor Uij = (diUj +djUi)/2 and we have used 
the elastic constants Cjj as defined in || . The equilibrium conformation will minimize the 
elastic free energy while maintaining a force-free configuration. 

First we consider the case just above H m where cqq vanishes, the flux-liquid @. In 
this case, when the flux-lines are aligned by an external magnetic field, we can directly 
show that a current parallel to the field tends to twist the flux-line tangents, as in a 
cholesteric liquid crystal. We employ the duality mapping between the superfluid and the 
superconductor under an applied local current. We write j° = p e ev° where v is the 



Cooper-pair velocity, p e is the pair density and e is the pair charge. The partition function 
for the London theory in an applied field H is: 

Z = J [dv][dA] <5[V- A] exp j- J d 3 x ^^(v - v° - A) 2 + i(V x A) 2 + H- V x A J , 

(10) 

where m is the mass of the Cooper-pair. Writing the velocity in Fourier space in terms of 
longitudinal and transverse components, 



/in , , ikxm 
v(k) = ik(t) + —— r 



'II) 



where m is the density of flux-vortices pointing in the m direction. Since flux-lines cannot 
begin or end in the sample V-m = 0. Upon substituting ([11]) into ( |T0| ) and integrating out 
4> and A we have (to leading order in momentum) 



[dm] <5 [V-m] exp 



drx 



1 



-m 



m-H- v°-V x H + v°-V x m 



(12) 



where we have omitted terms independent of m and v . The first two terms are responsible 
for the presence of vortices in the superconductor - they favor m = H. The next term 
induces the screening current in the Meissner phase. The last term is the new interaction 
which tends to twist the vortices around the applied current. If H is along z, it is natural 
to write m pz+po5n where 5n is the projection of the average tangent onto the xy-plane. 
Then V-m becomes: 

d z 8p + p V ± -5n = 0. (13) 

This constraint can be solved || by introducing a two-dimensional vector field u and 
writing 5p = —poS7^-u and 5n = d z u. In terms of this field u 



total 



(14) 



where we have allowed for anisotropic elastic constants. This theory is simply the theory 
of polymer cholesterics, in the limit of small pitch 0. For large deviations one might 
expect that the flux-lines will rotate in a plane perpendicular to a pitch direction. This 
configuration was, in fact, proposed in the seminal work of Campbell and Evetts 0. More- 
over, an additional conformation is possible in a finite radius, cylindrical sample, namely 
a double-twist configuration as in the blue-phase of chiral liquid crystals. This possible 
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double- twist configuration is shown in figure 1. In this double- twist conformation the B- 
field and j wrap around each other, simultaneously satisfying the Maxwell and London 
equations (|2|) and (P. It was correctly noted in [Q that this configuration would be ener- 
getically unacceptable as the radius of the cylinder grew - the flux lines on the boundary 
of the sample would grow unacceptably long. However, as we shall see in the following, a 
defect-riddled state can allow local configurations similar to those shown in figure 1, with 
finite displacements of the flux-lines. The handedness of the rotation of the flux-lines is 
determined by the right-hand-rule and the direction of the current. Note, however, that an 
equally acceptable force-free conformation would be absolutely straight flux-lines parallel 
to the applied current. The difficulty with this is that thermal fluctuations will destabilize 
this state and lead to a helical instability of flux-lines, as predicted by Clem [H in 1977. 
As a consistency check, we note that if we were to consider the effect of the Lorentz forces 



acting on the individual flux-lines that there is an instability also at any finite current [jlO 
towards helical flux-line trajectories. 

We note that there is a certain duality between the current and the magnetic field in 
the London- Maxwell equations. In particular, the equations are invariant under 

B -> j. 

c 

It would thus be natural to consider the dual physical situation to the Abrikosov flux- 
lattice. In this case, the current would flow along the z-axis leading to a screening magnetic 
circulation in the xy-plane. If the current were confined into regions, so would be the 
circulating magnetic field. Physically, this is accomplished via flux lines tracing out helical 
trajectories: the xy-components of the flux-line tangents circulate in that plane, dragging 
the magnetic field with them providing the necessary magnetic field. 

There is, however, an essential difference between the Abrikosov solution and its dual: 
in the original problem, quantum mechanics imposes a constraint on the amount of mag- 
netic flux that could be confined in a a flux-tube - single- valuedness of the wavefunction 
implies that the flux must be an integer multiple of the flux quantum $o = 2 2 ■ This 
constraint is responsible for the presence of a second-order transition between the Meiss- 
ner state and the Abrikosov state. In the dual case, there is no equivalent quantization 
of current flux. It is easy to understand why in the dual language: a helical flux-line can 



execute an arbitrarily long-pitched wobble which allows [V x B] i to be arbitrarily small. 
This is what allows Clem's helical instability. 

When / we are forced to consider a crystalline structure with a force-free 
conformation of the flux-lines. This problem has been considered in the context of liquid 
crystals - namely, how a chiral line-crystal minimizes its free-energy in the presence of 
the two-competing tendencies to twist and to have periodic order. These two tendencies 
frustrate each other and thus, as in the Renn-Lubensky twist-grain-boundary (TGB) phase 
of smectic liquid crystals the frustration will be resolved via the introduction of topological 
defects - screw dislocations. It is amusing that the TGB phase is the analog of the 
Abrikosov phase of the superconductor and so in the problem at hand we are minimizing 
the stresses of the real Abrikosov lattice with a "dual" Abrikosov lattice. 



In [|TTJ| two types of crystal defect arrays were considered. One array consisted of a 
periodic arrangement of tilt-grain-boundaries (TGB) which would change the local flux- 
line direction. The other array was made of helicoidal-grain-boundaries (HCB), each of 
which is a honeycomb lattice of screw-dislocations lying in the a;y-plane. A single isolated 
HCB leads to a twisting of the crystalline order along the flux-line direction. If we were 
to consider stacking many HCB's together with some spacing d', this twisted moire state 
would have both twisting of 5n as well as twisting of the crystal directions. This is similar 
to the physics of blue phases in chiral liquid crystals [|12|]. In chiral liquid crystals, there is 
a tendency for the local director n to twist. However, in blue phases this twist manifests 
itself in double-twist cylinders. Taking the nematic director field as a local tangent vector 
density for lines, these double-twist cylinders become rope-like bundles of twisted lines. 
Analogously, a twisted bundle of flux-lines will allow the magnetic field to circulate while 
keeping the flux-lines, on average along a single direction. While in the softer liquid crystal 
theory the elastic energy cost of this deformation is proportional to the angle of rotation 
PHI , in the flux-line system interactions between the screw dislocations of the vortex lattice 



will lead to logarithmic corrections to this energy. In any event, the energy of a grain 
boundary per unit area will be finite. We propose these defected states as paradigms for 
a flux-line lattice under an applied, parallel current. 

Notice that we can have no twist if a = 0. However, j = ctB. If current flows through 
the superconductor then it must flow on the oft-neglected boundaries of the sample. Thus, 
to study the energetics of this state, we must include the usual London energy for the 
supercurrent. We consider a current along the magnetic field direction z and assume that 



there is a flux-line lattice. If v is the Cooper-pair velocity and p c is the density of Cooper- 
pairs, then the total current is I z = jep c Av where A is the cross sectional area of the 
region in which current flows. If there are no defects in the flux lattice the current must 
flow within a penetration depth A of the sample boundary and A pa 2nRX. The London 
energy for this current configuration is, per unit length along z, 



R 



F/L = 2n I rdrm ePc v 2 = ™ e 2 I 2 Z , (16) 
Jr-\ 4TvR\p c e 2 

where m e is the mass of the electron. 

If we allow the Abrikosov lattice to have defects then the current can flow through 
more of the cross section thus lowering the London energy. Of course, the energy decrease 
will be offset by the energy of the screw dislocations in the flux-line lattice. If we consider 
a moire configuration which is reasonably dense then as a rough approximation we take 
V x Sn to be uniform along z . This implies that the current runs uniformly through the 
entire cross-section of the sample. The decreased London energy is: 

Fdefects/L = ^^II (17) 

In the moire state we must add the energy of the dislocation lattice that produces the 
twisted configuration. In a crystalline lattice the energy in the strain field due to a single 
dislocation diverges logarithmically with system size. If we have a network of dislocations 
as shown in figure 2, however, the strain energy of the lattice is finite. In this case the 
largest contribution to the energy of a grain-boundary is due to the logarithmic interactions 
of the screw-dislocations. We take the energy per unit length of a screw-dislocation to be 



eoln((i/C) where d is the average defect spacing and ( is the defect core size [fT3|j . The 
energy cost per unit length along z of the sample is therefore: 

7ri? 2 

^defects/ L = e ln(d/C)-^-. ( 18 ) 

With a uniform current density Maxwell's equation gives VxB = 4I z /(cR 2 ). In turn, the 
defect density is determined by the amount of V x 8n required to produce the uniform 



current. For a defect spacing d we estimate [14 



(19) 



Putting this together with Maxwell's equation, we get 

ttR 2 8irl z 1 



d 2 



(20) 



C y/$ H z 

where we have used po = H z /<& and ao = 1/y/po. Thus the energy cost from the defects 
is: 

E deiects /L = In V^f 2 ] 4 : e0h . (21) 

Putting together all the energies we may compare the energy of the moire structure 
with that of the untwisted Abrikosov structure. There will be an instability towards a 
twisted state when 



in 



-ii> 



in 



-J 2 + ^defects /L. 



(22) 



2nRXp c e 2 2 ~~ TcR 2 p c e 2 
If A <C R we may neglect the first term on the right hand side of fl2~2"|) and find that the 
moire state is favored when 



I z >ln 



Cy/$ H z nR 
SnI z C 



21 



8n 2 RXp c e 2 e 
m e cy/$ H z 



(23) 



or when the applied current density, j z = I z /(ttR 2 ) is 



A , 

Jz ~R ln 



Cy/$ H Z 

8nj z ( 2 



87rp c e 2 e 



m e cy/$ H z 



(24) 



Thus, as the system size increases the current density necessary to go to the moire state goes 
to zero. Thus the instability of Clem may be stabilized through the lattice structure and 
its screw dislocations. The moire state that proposed here is not exactly force-free. In each 
cell of the honeycomb, the flux-line displacement is a z-dependent rotation: m = CzeijXj. 
This configuration has V x B parallel to z, not the local B. However, when averaged over 
one cell, B is parallel to z and so there is no net force on each bundle of vortices. It is 
clear, however, that small adjustments to the flux-line locations that do not change their 
topology can yield an entirely force-free configuration. 

We end by commenting on recent experiments [|TS] performed on thin-film, 
YBa2Cu30y_5, a6-plane superconductors which have seen what are called "vortex"- 
twisters, composed of thousands of flux-lines. These twisters are formed by first applying a 
magnetic field along the long-direction of the sample creating an Abrikosov flux-line-lattice. 
An additional field H± is applied perpendicular to the plane of the flux-lines. As usual, 
the superconductor responds with a screening current which generates Lorentz forces on 



7 



the vortices. Because of the anisotropy of the superconductor the flux-lines will remain, for 
the most part, in the a6-plane and will not realign with the net magnetic field direction if 
H± is small enough |16| . On the short side of the sample these Lorentz forces push the flux 
lines along the direction either parallel or anti-parallel to the c-axis depending on the sign 
of the force. This bending will only occur near the sample edge as there is a competition 
for each flux-line between the in-plane pinning and the tendency to align with the net field 
\\17\\ . The forces build up and the Bean critical state []18[ is formed. For current to flow 
along the sides of the superconductor, the flux-lines must twist about, giving a non-zero 
V x B. Moreover, these vortex twisters can be made ever more stable and compact through 
"work- hardening" via an AC component of H±. When the AC field is removed, the work- 
hardened bundles are stable for hours, while the unhardened bundles are not stable at all. 
Since the moire state is a highly-ordered structure, it is unlikely that these vortex-twisters 
form a complete lattice of honeycomb-boundaries. However, it is possible that the local 
structure of these vortices resembles the highly entangled moire state. This could explain 
the work-hardening: when there are many screw dislocations it is difficult for the flux-line 
lattice to relax, as the defects cannot cross the flux-lines without cutting them. A flux 
bundle twisted as in figure 1 could relax easily as there are no topological impedements. 

Along these lines, it would be interesting to study (experimentally and theoretically) 
the dynamics of the flux-line lattice under an applied current. In particular, just above H m , 
the flux-line picture would suggest that there should be a striking frequency dependence of 
the I — V curve for currents along the field axis. If the frequency is smaller than a typical 
vortex diffusion time then the flux-lines can adjust their configurations to allow for current 
flow. However, as the frequency is increased the flux-lines will not be able to adjust for 
the current. This is reminiscent of the visco-elastic behavior of polymer melts. Moreover, 
if work-hardening of the Abrikosov lattice is possible, one might be able to increase the 



critical current in the sample ||19|| . These possibilities would verify the dual picture of 
high-Tc superconductors which focuses on the Abrikosov flux-lines. 
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sky, M.C. Marchetti, L. Radzihovsky and especially D.R. Nelson. We also thank the 
Aspen Center for Theoretical Physics, where some of this work was done. This work was 
supported, in part, by an award from Research Corporation and an NSF Career Award 
through Grant Number DMR97-32963. 
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Figure Captions 



Fig. 1. Configuration of flux-lines and current in a wire (both follow the heavy lines). 

The flux-lines are parallel to the current everywhere and both wrap around the 
center of the wire. Note that there is a nonvanishing V x B and V x j. The 
applied field H and the average current density j are parallel. 

Fig. 2. Proposed braided (moire) state of flux-lines. The dark lines make up the honey- 
comb network of screw dislocations. This texture has an approximately uniform 
V x B and so the current will flow uniformly through the cross-section of the 
sample. 
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Figure 1 



Figure 2 



